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Abstract 



We study the evolution of quantum eigenstates in the presence of level crossing under adiabatic cyclic change of environmental 
parameters. We find that exotic holonomies, indicated by exchange of the eigenstates after a single cyclic evolution, can arise from 
non-Abelian gauge potentials among non-degenerate levels. We illustrate our arguments with solvable two and three level models. 
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Berry phase phenomena are known to be wide-spread, in- 
triguing, and useful in controlling quantum systems 0~1|2). The 
Wilczek-Zee variation, in which different eigenstates sharing 
a degeneracy are turned into each other after the cyclic varia- 
tion of environmental parameter y), is found particularly use- 
ful, since it supplies the basis for so-called holonomic quantum 
computing [4|. Whether such state transformations require the 
existence of degeneracy throughout the parameter variation, is 
a matter in need of further analysis, although it has been the 
wide-spread assumption. 

The exotic holonomies are defined as exchange of eigen- 
states after one period of cyclic parametric evolution without 
any relevant degeneracy. They have been found in time-periodic 
systems |6) and in singular systems J7] [8], but not in fi- 
nite Hamiltonian systems up to now. One obvious reason is 
that the two real numbers in a real axis cannot be continuously 
exchanged without colliding with each other, i.e. degeneracy. 
Note that in a time periodic system described by a unitary ma- 
trix, its eigenvalues are complex number on the unit circle in the 
complex plane so that the eigenvalues smoothly exchange their 
position without crossing over each other. Once level cross- 
ing is allowed during the variation of environmental parameter, 
even in Hamiltonian system there is a possibility of the eigen- 
levels being exchanged after the cyclic parameter variation even 
among non-degenerate levels, which would be instrumental in 
enlarging the resource for holonomic quantum control. 

It appears that this is related to a myth that adiabatic param- 
eter variation excludes level crossing, and therefore, there is 
no possibility for exotic quantum holonomies for non-singular 
Hamiltonian system. An argument often raised against systems 



with level crossing is that it is not generic, and it represents a set 
of measure zero in parameter space of all systems. However, 
systems of our interest often lives in a world with some sym- 
metry, exact or approximate. Under certain circumstances, a 
system can always pass through the point of symmetry when an 
environmental parameter is varied along a circular path. Then, 
the level crossing become a common feature rather than an ex- 
ception. 

Surprisingly, it has been known for quite sometime that the 
adiabatic theorem is extendable to the case of crossing levels 
|9|. The eigenstates change smoothly as functions of environ- 
mental parameter even when the level crossing takes place, and 
this opens up the possibility for eigenvalue holonomy for non- 
singular Hamiltonian systems. In this article, we explicitly con- 
struct such models, which seem to have immediate extension to 
N level cases. The solvability of our model allows us to ex- 
amine the analytic structure of the gauge potentials, which is 
known to be the mathematical origin behind the existence of 
exotic holonomy IfTOlfTTII . 

2. Adiabatic level crossing and exotic holonomy 

Consider a Hamiltonian system with an environmental pa- 
rameter, which we call 6. We assume, for the moment, that the 
eigenvalues E n (ff) and eigenstates ^(fl) are non-degenerate for 
all possible values of the parameter. Let us suppose two levels 
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Figure 1 : Schematic diagram showing adiabatic level-crossing. A and D repre- 
sent energy gaps. 
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very closely approach each other at a certain value of the pa- 
rameter, which can be thought of as an avoided crossing with 
the closest energy gap A. Let us further assume that the two lev- 
els are separated apart from the next closest level (See Fig. 1) 
by D. Consider a smooth cyclic variation of 0(t) with a period 
t, namely, 0(t) - 6(0). Specifically, we require that 0(t) starts 
smoothly at the beginning t — and stops gently at t — t lfT2l to 
ensure the applicability of the Landau-Zener formula lfT3l[T4l . 
Let us assume that we have inequalities, 
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The period t is small enough compared with the inverse of the 
energy gap A so that the levels completely cross over the gap 
during the parametric variation, while it is large enough to ig- 
nore any transition among levels except the interacting two lev- 
els considered here, which is the reason why we call this pro- 
cess adiabatic. The situation remains intact even when the lev- 
els do cross according to some exact symmetry, A = 0, instead 
of showing avoided crossing. In fact, although a tiny avoided 
crossing caused by a slight symmetry-breaking takes place, the 
afore mentioned adiabatic level cross-over is robust irrespective 
of small parametric perturbation. This is the physics behind the 
so-called adiabatic level crossing fl9). 

Once such an adiabatic level crossing occurs, there is no 
reason to assume that each quantum eigenstate should come 
back to the corresponding initial state after a cyclic parametric 
variation. Only requirement is that the entire set of eigenstates 
should be the same as before, since the solutions of eigenvalue 
equation with a given parameter are uniquely determined. It 
is thus allowed that the two levels are exchanged after the para- 
metric variation. With a parameter 0, which describes the cyclic 
parametric variation along the path C from 0, to 6/ satisfying 
H(6j) - H(6f), such a transition is described by the holonomy 
matrix M as |[T5l[T6l 



Mfijn *¥ m (0i )e 



with 



M = T*e~ 



f c d9A(B) Te ij c d9A D (e) 



(2) 



(3) 



where T and T* represents the path-ordering and anti-path or- 
dering of operator integrals, the A{8) is the non-Abelian gauge 
potential 



and A D {6) its diagonal reduction 

A D nm {0) = A, un {0)5 n , m . 



(4) 



(5) 



The dynamical phase <p„ depends on the precise history of the 
parameter variation, while the holonomy matrix M is solely de- 
termined by the geometry of the path C in the parameter space. 
Non-zero off-diagonal component of M, if any, signifies the ex- 
istence of exotic holonomy. Physical requirement that an eigen- 
state does not split with adiabatic parameter variation limits the 
form of M to be permutation matrix supplemented by possible 



Manini-Pistolesi off-diagonal phases IfTTl fTBl for each non-zero 
elements. Namely, there is only a single non-zero entry to each 
raw and each column, and the absolute value of this entry is 
one. 

3. Two-level model with exotic holonomy 

Consider a two level quantum system described by a para- 
metric Hamiltonian 



H(0) = R(ff) 
with a real number v, 

Z® =a 7 = 



cos -Z 
2 



(2) 



vsin —F 
2 



(2) 



1 
-1 



F w = p +(rx = 



and anti-periodic function R(0) with period 2n, 
R(8 + 2n) = -R(0). 



(6) 



(7) 



(8) 



A convenient choice we adopt in the numerical examples is 
R{0) = cos |. The system thus becomes 2n periodic; 

H{0 + 2n) = H(8), (9) 

and the parameter e [0, 2n) forms a ring, S 1 . Note that if ^ 
is written in the form H(x,y) = xZ <T> + y vf (2) , the parametric 
evolution : — > 2n represents a circle on (x, y) plane with its 
center shifted by the radius into x axis, namely x — 1 + cos 8 
and y = sin 0, so that it touches the origin. If we were to vary 
R e (-00,00) and e [0,2^) independently, the entire plane 
(x, y) is covered. 

The eigenvalue equation 



HQEWM = E n (0f¥n(e) 
is analytically solvable with eigenvalues given by 



E„(0) =R(6) cos -P„(6), 



and eigenstates 



P„{0) + 1 



y/2P„(0) 2 + 2 ~ 1 

for n — 1,2, with 



(10) 



(ID 



(12) 



P„(0) = vtan-+(-)"sgn 



V cos ■ 



1 +v 2 tan 2 -. (13) 



Energy eigenvalues as a function of are shown in Fig. 2, in 
which the most notable feature is the occurrence of degeneracy 
at — 7i and the related exchange of eigenvalues. It guarantees 
that the set of eigenvalues at = is equivalent to that at — 2n, 
which is a direct result of H(0) being identical to H{2n). The 
degeneracy of eigenvalues at 9 — n is a direct consequence of 
vanishing Hamiltonian, H(n) = 0. 

We note that P„(6) is 47r-periodic. Moreover, we have Pi (0+ 
2tt) = P 2 (0) and P 2 (0 + 2n) = P\(0), which leads to the appear- 
ance of exotic holonomy, i.e., 



Ei(0 + 2tt) = E 2 (0), E 2 (0 + 2n) = E X {B), 



(14) 
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Figure 2: Energy eigenvalues E„(6) of the model (6) as a function of 8 with 
J_ 



v 



where 



1 d X (ff) 
4 39 



(19) 



In Fig. 4, we depicts two examples of function f(9). Obviously, 
we have C d9f(9) = | , and we obtain the holonomy matrix 



M : 



1 
-1 



(20) 



showing the exotic holonomy with Manini-Pistolesi phase (- 1 )" . 



and also, *Fi(6> + 2n) oc ¥ 2 (0) and ¥ 2 (0 + 2tt) oc ¥,(0). 

The structure of the eigenstates becomes clearer with the 
re-parameterization of P„(6) with new angle variable^ = ;^(0), 
which we define as 



P n (X) = tan 



X + (2n - 3)n 



(15) 



namely, P\(x) = tan and f^Ot) = tan . The monotonously 
increasing function ^ (0) maps e [0, 27r) to ^ € [0, 27r). The 
eigenstates is written, with the new angle parameter^, as 

/ . x -(2n-6)n \ 



namely 

T,(Y) 



- cos 



COS ' 



sin 



(17) 



Note that ^,,(0) is htn-periodic, with anti -periodicity of period 
47T. This is not immediately evident from the expression (12) 
which is in fact discontinuous at several values of — nn, and 
had to be amended with factor -sign(sin " ) to turn into 



(16 1. A numerical example of eigenstates as functions of is 
depicted in Fig. 3. 




Figure 3: Eigenstates (bottom) and ¥9(0) (top) of the model (6) with 

R(8) = cos j , and v = ^4= . The solid and the dashed lines represent the up- 
per and the lower component of the eigenvectors, respectively, both of which 
are chosen to be real. The range indicated by thick lines represents a single 
period e [0, 2n] . The values outside of this range are shown to display the 
periodicities and mutual relations of *P|(0) and I^Cff) 

Adiabatic change of eigenstates are determined by the gauge 
potential A nm (8) given by 



A{6) = 



(18) 
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Figure 4: Functional form of gauge potential f(6) of the model (6) with v 
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Figure 5: Exceptional points on the Mercator projection of the complex 9 plane 
of system described by (6) with v > 1 (left) and v < 1 (right). The filled crosses 
represent the exceptional points that are the poles of gauge potential A(0), while 
the unfilled crosses are the points of eigenvalue degeneracy which has no effect 
on the singular behavior of A(6). The solid lines are the branch cuts on which 
Re(ii2 - Ei) = 0. In the limit v — > 1, the two complex exceptional points (the 
filled crosses) move to ±ioo. 



Nontrivial holonomy is known to be related to the analytic 
structure of the gauge potential A(6) in the complex 9 plain 1 10 1, 
specifically, its singularities. These singularities can arise at the 
exceptional points 9* , which is defined as the point where two 
complex energy coincide; E m {9*)-E n {9*) = 0. In our example, 



R(9*) cos — (P 2 (9* ) - = 

which yields 

9* = n , 9* = 2arccot (+;' v) , 



(21) 



(22) 



where the first one coming from R(9*) cos y = 0, and the 9* 
from P2(9*) - Pi(9*) = 0. The exceptional points coming from 
R(9*) cos y = has obviously no bearing on the singularity of 
A(9), while at 9* we have the poles of A(9) in the form 



A 12(0) = -A 2l (0) * + 



2 1 



49-9* 



(9 -» 9* ± ). 



(23) 



The existence of the poles explains the non-vanishing values 
of j>d9A{2(9) and j>d9A\2(9) around the real axis, i.e. M12, 
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M%\ + implying the existence of exotic holonomy. Fig. 5 
shows the locations of poles and branching structure of energy 
surface E„(9) in Mercator representation of complex parameter 
plane 9. 

4. Three-level model with exotic holonomy 

Let us now consider a three level quantum system described 
by a parametric Hamiltonian 



H(9) = R(9) 
where v is real, 



cos -Z (3) + vsm-F (3> 
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and R(9) is anti-periodic with period 2jt, i.e. 
R(0 + 2n) = -R(9). 



(24) 



(25) 



(26) 



As before, in the numerical examples, we adopt R(9) - cos | . 
The system then becomes 2n periodic; 



H(9 + 2n) = H(9), 
and the parameter 9 e [0, 2n) forms a ring, S 1 . 



(27) 



E o 




Figure 6: Energy eigenvalues E„(8) of the model (24) with v = 1. Exotic 
eigenvalue holonomy is clearly observed. 



The eigenvalue equation 

H(ey¥M = E n (9w„(0) 

is analytically solvable with eigenvalues given by 



E n (6) =R(d) cos -Q n (6), 



and eigenstates 



1 



V3G«(0) 4 + 1 

for n - 1, 2, 3, with 

sgn(cos^) 



(Qn(8)(Qn(ff) + D) 

Q,M 2 - 1 

Qn(0)(Qn(O) - 1)J 



(28) 



(29) 



(30) 



Qn{0) = 



^3(1 - sini 2f ) 



. i rffl . rj{ff) - 2rj n 

sin 3 2 sin 

2 6 



(31) 



in which the angle rj{9) is defined by 
j](9) = 2 arcsin 



3v 2 tan 2 f 



1 + 3v 2 tan 2 



2 > 



(32) 



and the state dependent shift r\ n = (2n - 4)n. The function rj(9) 
is monotonously increasing and maps 9 e [0, 2n) to 77 e [0,2^). 
One example of the energy eigenvalues as function of environ- 
mental parameter 9 is shown in Fig. 6. All eigenvalues are de- 
generate at 9 — 7i, as a consequence of vanishing Hamiltonian 
H(tt) = 0. 

We note that Q„(9) is 67r-periodic. Moreover, we have Q„(9+ 
2n) - Q n +\{9), where the subscripts are to be understood in the 
sense of modulo three, which clearly signifies the existence of 
exotic holonomy, 



E n (9 + 2n) = E n+ i(9), (n = 1,2,3). 



(33) 
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Figure 7: Eigenstates (bottom), W 2 (0) (middle) and *P 3 (6) (top) of the 

model (24) with R(0) — cos and v — 1. The solid, the dotted, and the 
dashed lines represent the upper, the middle, and the bottom components of 
the eigenvectors, respectively, all of which are chosen to be real. See also the 
caption of Fig. 3. 

The structure of the eigenstates becomes clearer with the 
re-parameterization of Q„(9) with a new angle variable £ = £(9) 



Qm = Rt 



(34) 



namely, = Rt [(f - 2n)/6], Q 2 (0 = Rt [f/6] and Q 3 <£) = 

Rt [(£ + 2n)/6], where a ^-periodic function Rt is defined by 

Rt f = — sgn(tan Vltanfl = , (35) 
33 Rc£ 

along with 27r-periodic functions Rs and Rc defined as 

Rs£ = —sgn(sin^) Vising, 
3? 

Rc £ = sgn(cos yj cos £|. (36) 

The functions Rt, Rs and Rc are analytic on a single sheet 
complex 9 plane in contrast to Vtan -y/sin and a/cos£, re- 
spectively, which are analytic on a double sheet 9 plane. The 
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monotonously increasing function £(9) maps 9 e [0, 2n) to f € 
[0, 2n). The eigenstates is written, with the new angle parame- 
ter as 



Rs^(Rs^-Rc^ 



Rs 



2 f+Ift 



„ Rc 2 ^ 

Rs ^ (Rs ^ - Rc ^ 



)J 



(37) 



From this from, we see that v i'„(9) is 6n-periodic. 

The gauge potential A„ m (6), which determines the adiabatic 
variation of eigenstates, is given by 



(38) 
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*(0-2w), 



where g(ff) is defined by 
6 30 



< ( P 1 (^))|/3^3(^))>. 



(39) 



The calculation of the holonomy matrix involves fully ordered 
matrix integral, thus no simple analytical calculation can be per- 
formed. However, we can deduce from ((37), that it is given by 



M = 



1 

1 

1 



(40) 



showing the spiral type exotic holonomy with Manini-Pistolesi 
phase 1 for all states. 




Figure 8: Functional form of gauge potential g(6) (solid line), g(fi— 2n) (dashed 
line) and g{8 + 2ri) (dotted line) of the model (24) with v = 1 . 

As in the two level case, we examine exceptional points 9* 
where two complex energies coalesce, E m (9*) - E„(9*) - 0. 
We find 



(12) 



n*(23) 



(41) 



where the first solution comes from cos y = 0, while f?* (12) 
and 0* (23) are obtained from P 2 (0*) - Pi(0*) = and P 3 (9*) - 
Pi(6*) = 0, respectively. We then immediately obtain 



k(12) 



♦•(23) 



-2 arccot \3(e 



2 arccot 



^(e***- l)v 



(42) 



Near the exceptional points, A(0)s are approximated as 

i 1 



A Jk (6) = -A jk (8) 



(9 -> <9* w ). (43) 



The existence of the poles explain the non-vanishing values of 
M\z, M23, and M31, implying the existence of exotic holonomy. 
Fig. 9 shows the locations of poles and branching structure of 
energy surface E„(9) in Mercator representation of the complex 
parameter plane 9. 
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Figure 9: Exceptional points on the Mercator projection of the complex 8 plane 
of system described by (24). The filled crosses represent the exceptional points 
that are the poles of gauge potential A(8), while the unfilled cross is the point of 
eigenvalue degeneracy having no effect on the singular behavior of A(8). The 
solid and the dashed lines represent the branch cuts satisfying Re(£3 - E2) = 
and Re(i?2 - E\ ) = 0, respectively. 

All the results obtained here do not depend on the specific 



choice of the matrices (25 I. In fact, we can enlarge our model 
by replacing the two matrices Z (3) and by 

Z (3) =S4 + c 5 E5, ^ (3) = 7 (3) +C1Z1 +C2E2 +c 3 23, (44) 

where cj are real numbers, 7 (3) three-dimensional unit matrix 
and Ej given by 
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(45) 



If we were to make independent choice of six parameters, R e 
(-00,00), 9 € [0,2n), cj e (-00,00), (j = 1,3,4,5), the Hamil- 



tonian ( 24 1 is the most general real symmetric three-by-three 
matrix. By fixing c,s and binding R and 9 by R{9) = cos I, 
we go back to the same game of considering the system as a 
function of a single parameter 9 which forms a ring S 1 . It can 
be checked numerically, that the exotic holonomy characterized 
by (40i, or equivalently, the eigenvalue flow {1, 2, 3} — > {2,3,1}, 
in obvious notation, is a common characteristics of system with 
c\ — C2 — C3. Here, we have made an assumption that the unper- 
turbed spectrum is not much different from the original model, 
|<75 1 <sc 1, All possible patterns of eigenvalue flow are obtained 
with suitable choice of cjS. Specifically, c\ + 0, C2 — C3 = 
results in {1,2,3} -> {2, 1,3), c 2 + 0, c 3 = c\ = in {1,2,3} -> 
{3,2,1}, and c 3 * 0, c x = c 2 = 0, in {1,2,3} {1,3,2}. 
Generic case c\ + c 2 + c 3 also produces the second pattern, 
{1, 2, 3} — > {3, 2, 1). It is now clear, that there is finite subset of 
parameter space, in which exotic holonomies of various types 
arise after cyclic variation of the parameter 9. 



5. Outlook 

Our results obtained in the two and three level cases can 
be extended to N levels (N < 4) in a straightforward way. It 
is possible to prove the existence of the exotic holonomy for 
systems described by the Hamiltonian 



H(8) = cos 



cos -Z' 
2 



. 
+ v sin -\w 



(N) 



)(w' 



(,<V)i 



(46) 



where |w^) is a normalized A^-dimensional vector and Z (AP) is 
an N x N Hermitian matrix, as long as all eigenvectors of Z (AP) 
have non-zero overlap with |w (Ar) ). 

The Hamiltonian exotic holonomy shares a common feature 
with the Wilczek-Zee holonomy of having SU(N) non-Abelian 
gauge potential at their base. However, they are distinct in that, 
in the former, N eigenstates are exchanged among themselves 
with their internal dynamics, while, in the latter, involvement 
of another eigenstate, or a set of degenerate eigenstates j3) is 
required. 

The exotic holonomy also seems to have resemblance to 
the off-diagonal holonomy of Manini and Pistolesi. It is im- 
portant to point out that the eigenstates are obtained indepen- 
dently from the choice of envelope function R(6). If we make 
the choice R(ff) - 1, the Hamiltonian becomes anti-periodic, 
H(8 + 2n) = -H(8), so that new period is now An. The eigen- 
state holonomy, occurring now at the midpoint of the new full 
cycle 8 e [0, An), is nothing but the off-diagonal holonomy dis- 
cussed by Manini and Pistolesi. In the off-diagonal holonomy, 
the set of eigenstates at the starting value of environmental pa- 
rameter "accidentally" coincides with that at another value of 
parameter in a midpoint of cyclic evolution. In general, such 
coincidence is highly unlikely, and it is often a result of the 
same Hamiltonian multiplied by different numbers appearing at 
different value of environmental parameter. In such a case, with 
the introduction of new envelope function and reinterpretation 
of the period of parameter variation, a system with off-diagonal 
holonomy can be mapped to another one with exotic holonomy. 
In this instance, the Manini-Pistolesi holonomy is the exotic 
holonomy in disguise. 




parameter space, whose connected apices of two cones repre- 
sent Berry's diabolical point. When the parameters are varied 
along a circle that surrounds the diabolical point, Berry phase 
arises (Fig. 10 left). We can ask a question: what will hap- 
pen when the circle touches the diabolical point. Obviously, 
the trajectory on the energy surface should be "smooth", and 
it wanders both cones. With a cyclic variation of parameters, 
the trajectory moves from one cone to the other. This is exactly 



the Hamiltonian exotic holonomy (Fig. 10 right). This process 
is fully described by holonomy matrix given in terms of the 
path-ordered integral of the gauge potential just as in the case 
of Berry phase. The gauge potential now has singularities in 
complexified parameter space, not on the diabolical point itself. 

The Hamiltonian exotic holonomy can be viewed as an ex- 
tension of, and a natural complement to the Berry phase, and it 
forms an integral part of physics of adiabatic quantum control. 
The general equation for quantum holonomy ([3]l is just a very 
natural expression of the basic requirement that the entire set of 
eigenstates is to be mapped to itself after the cyclic variation of 
environmental parameter. 
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Figure 10: Two types of Hamiltonian holonomy in the presence of a diabolical 
point (cross surrounded by small circle) on energy surface standing on a para- 
metric plane. The picture on the left represents circular parameter variation that 
results in the Berry and Wilczek-Zee holonomies, while the picture on the right 
depicts the one leading to the exotic holonomy. 



Our findings can be placed in context by considering a stan- 
dard double-cone structure of energy surface standing on the 
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